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Abstract 

A unified form of mass matrix proposed previously for all fundamental fermions is 
extended to include the sterile neutrino u s , presumed to mix mainly with u e , while v T mixes 
strongly with (and both and v T weakly with v e and v s ). It turns out that the former 
mixing can be responsible for oscillations of solar neutrinos, while the latter for those of 
atmospheric neutrinos. It is mentioned in a footnote that v s is one of two sterile neutrinos, 
possible from a viewpoint of an algebraic construction. The charged-lepton factor 
in the leptonic four-dimensional Cabibbo — Kobayashi — Maskawa matrix V = U^^U^ 
is determined through the small deviation of our previous prediction m T = 1776.80 MeV 
(valid at the level of = 1) from the experimental value m T = 1777. OOlg;^ MeV. Then, 
it gives corrections to the four-dimensional mixing of z/ e , z/ M , u T , v s leading, in particular, 
to small v e — > u^, v e — > v T and — > v s oscillations (absent at the level of = 1). 
However, the LSND estimates of — > v e and z/ M — > v e oscillations are too large to be 
explained by these corrections, though in the case of — > v e they seem to be at the edge 
of such possible explanation. 
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1. Introduction 



Recently, in the course of extended studies on the "texture" of charged leptons e~ , //~ , 
T~, neutrinos u e , u T , up quarks u , c, t and down quarks d , s , 6, we came to a proposal 
of unified algebraic structure of their mass matrices (M-j^ (/ = v , e , u , d) in the three- 
dimensional family space = 1,2,3) [1]. The proposed structure followed from two 
sources. First of all, from (i) an idea about the origin of three fundamental-fermion 
families as a consequence of some generalized Dirac-type equations (interacting with the 
Standard-Model gauge bosons) whose a priori infinite series is, due to an intrinsic Pauli 
principle, reduced (in the case of fermions) to three such equations^ [2]. And further, 
from (ii) an ansatz for the fermion mass matrix expressed in terms of the suggested 
family characteristics [2]. 

This proposal in the case of leptons reads 
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with f = v , e, while /i^ , e^ 2 , , and denote real constants to be determined 
from the present and future experimental data for lepton masses and mixing parameters 
, «W and pW are mass-dimensional). 
On the base of a numerical experience we assumed that 

e^ 2 = 0, a {u) = 0, /5 (e) = 0, (2) 

what leads to Mq = and = = M21 . Then, we determined the parameters //^ , 

e^ 2 and a* 6 - 1 from the experimental values of m e , and m T , while and (3^ - 

tThese imply the existence of three Dirac bispinors ipa , 4>a — (1/4)(C _1 7 5 ) / 3 1 /3 2 ip a 1 p 2 and ipa — 
(l/24)e / 3 l/ 3 2/ 3 3 / 3 4 ^ Q/ 3 lA/ 3 3/ 3 4 ,wherea= 1, 2, 3, 4 and & = 1,2,3,4 (z = 1, 2, . . . , JV-1 with iV-1 = 0,2,4) 
are Dirac bispinor indices (a is correlated with the presence of a whole set of Standard-Model charges 
suppressed in our notation, while f3\, . . . , Pn-i are fully antisymmetric). We interpret these bispinors as 
fundamental fermions from three families, corresponding to a given Standard-Model signature, i.e., as 
v e , Vp, v r or e~ , fx" , r _ or u , c , t or d , s , b. 
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from the possible atmospheric-neutrino oscillations as seen by Super-Kamiokande. The 
phases (p^ and ip^ remained undetermined. 

In fact, in the case of charged leptons, assuming that the off-diagonal elements of 
the mass matrix (m/J^ can be treated as a small perturbation of its diagonal terms, we 
calculate in the lowest (quadratic) perturbative order [1] 
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When the experimental m e and m M [3] are used as inputs, Eqs. (3) give [1] 
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We can see that the predicted value of m T agrees very well with its experimental figure 
m cx P _ i777.oo^o;2? MeV [3], even in the zero-order perturbative calculation. In order to 
estimate , we can take this experimental figure as another input. Then, 
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so it is consistent with zero. 

For the unitary matrix (jj-f^, diagonalizing the mass matrix (^M-f^j according to the 
relation = diag(m e , m M , m T ), we obtain in the lowest (linear) perturbative 

order in / /x^ 
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where the small e^ 2 is neglected. 

The case of neutrinos is discussed throughout the next Sections. In particular, we will 
determine the parameters and (3^ v > in Section 4 from Super-Kamiokande results. 

We shall not discuss any a priori motivation for the proposal (1), considering it simply 
as a detailed conjecture (the interested Reader may look for its roots in Refs. [2]). Instead, 
we allow in this paper for the existence of a sterile neutrino v s (blind to the Standard- 
Model interactions), extending the neutrino mass matrix [m\^ (i,j = 1,2,3), as given 
in Eq. (1), to a 4 x 4 Hermitian matrix (M^A (I, J = 1,2,3,4), where M$ = M$*f\ 
To this end, we supplement the former matrix by seven a priori unknown elements , 
M$ and Afff . We will assume by analogy with M$ = = and M$ = = M$ 
that 

M<? = = iwff , M<? = = M%> , (7) 

but allow for nonzero and M$ (as well as and Mj$) in analogy with nonzero 
M 2 ( 3 ] and M$ (as well as M$ and M$). Note that here = under our particular 
assumption (2). If there is no sterile neutrino v s , then also My± = = M$ and — 0, 
and we return to the 3x3 mass matrix {mI^ (i,j = 1, 2, 3) discussed in Ref. [1]. 

2. Neutrino mass states 

The eigenvalues of the mass matrix ^M^J^j are the masses of four neutrino mass states 
v \i v 2i ^3 , ^4- They are given as follows 
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* Roots for this sterile neutrino may be sought again in the generalized Dirac equations (this time, 
without Standard-Model interactions) whose a priori infinite number is, due to the intrinsic Pauli 
principle, reduced (in the case of fermions) to two such equations for the Dirac bispinors ipjjp and 
ipjp = (l/6)(7 5 C)p p i ep l 01 fa foil) fa fa fa {Pi , 02, fa are fully antisymmetric). In the present paper, ipjp 

(2) 

is interpreted as a sterile neutrino of the Dirac type denoted by v s . The existence of ipp as another sterile 
neutrino of the Dirac type is not discussed here (such second potential sterile neutrino might mix mainly 
with i/fi, producing an extra disappearance mode of v^, slightly correcting the effect of its dominating 
mode — > v T \ for some comments cf. Section 6). 
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The corresponding unitary matrix (Ujj^j, diagonalizing the mass matrix (M\ U j^j ac- 
cording to the equality = diag(m^ , m U2 , m Vz , m VA ) , takes the form 
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and |M 2 3 | 7^ for (3^ > 0. If there is no sterile neutrino, then Y — > as seen from Eq. 
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The neutrino states v a = u e , , v T , u s (of which z/ e , 1/^ , v T denote the familiar ob- 
served neutrino weak-interaction states, while v s stands for their unobservable sterile part- 
ner) are related to neutrino mass states vj = v\ , v 2 , v 3 , z/ 4 through the four-dimensional 
unitary transformation 
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with (Vj a ) = (F Qj ) t . Here, 
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with the charged-lepton diagonalizing matrix (Uif^j as given in Eq. (6) and 
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The latter equations follow from the fact that charged leptons get no sterile partner. 
Thus, 



Vaj = E U kT U kj , V aA = US2* . (14) 

k 

Of course, the 4x4 unitary matrix (V a j) is the four-dimensional lepton counterpart 
of the familiar Cabibbo — Kobayashi — Maskawa matrix for quarks. The charged leptons 
e~ , , t~ (with diagonalized mass matrix) are here counterparts of the up quarks u , c ,t 
(with diagonalized mass matrix). 

From Eqs. (12) as well as (9) and (6) we can calculate the matrix elements V a j in 
the lowest (quadratic) perturbative order in / (j,( e \ The result reads (we write for 
convenience a — I — 1, 2, 3, 4): 
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3. Neutrino oscillations 



Once knowing the elements (15) of the lepton Cabibbo — Kobayashi — Maskawa matrix, 
we are able to calculate the probabilities of neutrino oscillations u a — > vp (in the vacuum), 
making use of the familiar formulae: 



P(v n 



; m l L - m l K _t 



2IPI , 

\v a and (vp\v a ) = 8p a . Here, as usual, t/\p\ = L/E (c = 1 



(16) 



where v a (0) = v a , (u a \ = 
ft) and this is equal to 4 x 1.2QQ3L/E if ml — m 2 UK , L and E are measured in eV 2 , m 
and MeV, respectively. Of course, L is the source-detector distance (the baseline). In 
the following, it will be convenient to denote xlk = 1.2663(m^. — m 2 UK )L/E and use the 
identity cos 2xlk = 1 — 2 sin 2 xlk- 

From Eqs. (16) and (15) we deduce by explicit calculations the following neutrino 
oscillation formulae valid in the lowest (quadratic) perturbative order in a^/pL^: 
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for the appearance experiments (with the appearance modes of v T \ and 
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P {Vp -> = 1 - P {Vy, ^V e )~P (Vp -> I/ T ) 

16 /«( £ )\ 2 Y 2 2 

~84i (TTW sin * 41 ' (21) 

P(z/ r ^z/ r ) = l-P^^^-P^^^) (22) 

for the disappearance experiments (with disappearance modes of u e , u^, u T ). Of all these 
formulae, only Eqs. (18) and (19) have two-family forms. 

Notice that Eqs. (17) — (19) are invariant under the simultaneous substitution (pM — > 
—(pW and ip^ — > — <p( e \ what means their invariance under the replacement Vko. — > Vka- 
This implies through Eq. (16) that P(vp — > u a ) = P{y a — > vp). Thus, our neutrino os- 
cillation formulae, valid in the lowest (quadratic) perturbative order in / preserve 
T reversal and, by CPT invariance, also CP reflection, though the matrix (V a j) is here 
not real (the effect of nonreal (V a j) appears in higher perturbative orders and then spoils 
CP conservation). Note that CP conservation, when it works, and CPT invariance imply 
P(u a -> up) = P(u a — >• %) and P(u a -> vp) = P(Vp -> i> a ), respectively. 

Obviously, Eqs. (20)— (22) tell us that 
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for the "appearance" modes of unobservable sterile neutrino v s . 

If there is no sterile neutrino v s , then Y = wherever it appears in Eqs. (17) — (24). 
On the other hand, if (cc^/m^) 2 — > 0, then 
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4Y 2 

P {y e - „ c ) - 1 - (1 + y2)2 sin 2 X4i (28) 

and 

4X 2 

i 5 K ^ ^) ^ 1 " (1+X 2)2 Sin2 ^2 , (29) 

but P [y T — > i/ r ) — > 1. 

Concluding, we can see that the practically decoupled oscillations z/ M — > z/ T and v e — > z/ s 
play in the framework of our neutrino "texture" an exceptional, dominating role. If there 
is no sterile neutrino v s i.e., Y — 0, then such a role is played only by — > z/ r . 

Note that in Eqs. (17) — (24) sin 2 ^^ — > 1/2 if x LK /n — > oo, since the source and 
detector have always finite extensions over which sin 2 xlk ought to be averaged (with 
respect to their distance L). If 2xlk/' k 0, then sin 2 ^^ =— > x\ K — > 0. 

4. Information from atmospheric neutrinos 

The atmospheric neutrino experiments seem to indicate that there is a deficit of atmo- 
spheric i^'s, caused by the neutrino oscillations corresponding to disappearance modes of 
v^. These result in the survival probability which, if analized in two-family form 

P (^ -> Vv ) = 1 - sin 2 26 atm sin 2 (l.27Am 2 tm L/e) , (30) 

leads to 

sin 2 2# atm = 0(1) ~ 0.8 to 1 (31) 

and 

Am 2 tm ~ (0.03 to 1) x 10" 2 eV 2 (32) 

with the preferable value Am 2 tm ~ 0.5 x 10~ 2 eV 2 [4,5]. It is usually suggested that, 
practically, the oscillations — > u T alone are responsible for such a deficit. Then, 

P W -> v T ) = sin 2 20 atm sin 2 (l.27Am 2 tm L/P) . (33) 
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The last suggestion, neglecting the disapearance mode — > z/ e , is consistent with the 
negative result of CHOOZ long-baseline reactor experiment that finds no evidence for 
neutrino oscillations corresponding to the disappearance modes of z/ e , in particular v e — > 
[6]. The region of sin 2 2# atm and Am 2 tm indicated by Super-Kamiokande atmospheric- 
neutrino experiment [Eqs. (31) and (32)] lies, in fact, inside the region of sin 2 2#cH and 
Am^jj excluded by CHOOZ (where P (z7 e -> p e ) = 1). 

This important message from CHOOZ experiment, restricting the strength of mixing 

with v e , leaves a priori three options for mixing with v T and v s : mixes dominantly 
with v T (while v e with is s ), or with v s (while v e with is T ), or with both v T and v s (while v e 
does not mix). Of these three options, our neutrino "texture" chooses the first due to the 
assumptions = and M$ = 0, supplemented by M^ 1 = and m!$ = 0. Note that 
the second option would correspond to different assumptions, = and = 0, 

supplemented by M$ = and M$ = 0; eventually, in the case of the third option there 
would be M$ = 0, supplemented by M$ = 0, = and Afff = 0. 

We can see from the experimental atmospheric-neutrino estimates (30) and (33) that 
they correspond exactly to our neutrino oscillation formulae (29) and (26), respectively. 
Hence, we can infer that 



and so, with the Super-Kamiokande figures (31) and (32), we can take 
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0.8 to 1 



and 




5 x 10~ 3 eV 2 
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as two neutrino inputs. 

From the input (35) we evaluate the limits 



0.618 to 1. 



(37) 
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However, it is not difficult to see that (because the difference (36) is kept fixed) the 
limit of X — > 1 is singular in the sense that then ^ — > and f3^ — > oo with 
^M/JM _> (5/20.9) x 10- 3 eV 2 (and so, -> 0, Afff -> and |m£ } | -> oo with 

2(Mg } + M^M^I -> 5 x 10~ 3 eV 2 ). Then, =-> ^0.478 (3 {u) -> ^°o. We will 

restrict, therefore, the range in the input (35) to 

4X 2 

-0.8 to 1-10" 6 , (38) 



(1 + X 2 ) 2 

where the particular upper limit is chosen as an illustration. Hence, 

X~ 0.618 to 0.999. (39) 
In such a case, the second Eq. (10) leads to 



-(0.500 to 0.001) (40) 



M% - M&> X 2 - 1 



2\M&\ 2X 



'23 



(respectively). On the other hand, the mass matrix (m^^ as given in Eq. (1) implies 
that 

—22 33_ = _ 2 0.3^ rT . 41 

2\M^\ P {u) 
Thus, from Eqs. (40) and (41) we get 

^4 ~ 0.0246 to 0.0000492 , ^ ~ 40.7 to 20300 (42) 

(respectively). 

Then, making use of Eq. (1) for the mass matrix (M^^, we obtain from the second 
Eq. (8) 



m 



10.9^0.478-^. 1 + 

-(10.8 to 9700)/^) 
(32.7 to 9730)/i^ 



20.3 



/i 



(43) 



(respectively). Therefore, 
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ml 3 - m 2 = (951 to 425000)/i M2 



(44) 



Hence, using the input (36), we evaluate 

/i H ~ (0.00229 to 0.000108) eV (45) 
(respectively). With the values (45) of fj,^\ Eqs. (43) lead to 

m U2 (0.0247 to 1.05) eV , m U3 ~ (0.0749 to 1.05) eV (46) 

(respectively). Here, the minus sign at m U2 is phenomenologically irrelevant in relativistic 
dynamics (c/. Dirac equation). Note that 

m 2 V2 ~ (0.000611 to 1.11) eV 2 , m 2 3 ~ (0.0561 to 1.11) eV 2 , (47) 
where m 2 3 — m 2 2 ~ 0.005 eV 2 . Finally, from Eqs. (42) and (45) we evaluate 

(3 {u) ~ (0.0933 to 2.20) eV . (48) 

To summarize the above discussion of atmospheric neutrinos, we can say that the 
Super-Kamiokande experiment seems to transmit to us an important message about 
strong mixing of and v T and their rather weak mixing with v e . Such a situation is pre- 
dicted just in the case of our neutrino "texture" with X = 0(1) and small (ofo / — > 
0, the latter value being motivated by our excellent zero-order prediction of m T . This 
conclusion is independent of whether the sterile neutrino v s exists or does not exist in 
our neutrino "texture". The predicted neutrino masses are \m V2 \ ~ (0.02 to 1) eV and 
m U3 ~ (0.07 to 1) eV with m 2 3 - m 2 2 ~ 5 x 10" 3 eV 2 , where the upper limit X = 0.999 
is considered for X < 1. For X still nearer to 1, m U2 and m U3 increase further. 

5. Information from solar neutrinos 

As is well known, the solar neutrino experiments demonstrate a deficit of solar z/ e 's 
reaching the Earth, that seems to be caused by neutrino oscillations corresponding to 
disappearance modes of u e . These modes result in the survival probability, usually analized 
in two-family form 
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P {y e -> u e ) = 1 - sin 2 2# so i sin 2 (l.27Am 2 ol L/E) . (49) 

Here, the oscillation amplitude sin 2 26 so \ is likely to be enhanced to a new (sin 2 2# so i) mat t C r 
by the resonant MSW mechanism [7] in the Sun matter (dependent on values of sin 2 2# sol 
and Am 2 ol ), though the vacuum mechanism is still not excluded. According to recent 
estimations [8], in the first case there are two solutions with 

sin 2 2# sol ~8x 1(T 3 , Am s 2 ol ~5x KT 6 eV 2 (50) 

and 

sin 2 2# so i ~ 0.6 , Am s 2 ol ~ 1.6 x 10~ 5 eV 2 , (51) 
while in the second case 

sin 2 29 sol ~ 0.65 - 1 , Am 2 ol ~ (5 - 8) x 10" 11 eV 2 . (52) 

If the disappearance mode v e — > u s dominates, then from the two MSW solutions only the 
first survives. From the above three solutions, the first is considered as most favorable. 

Note that the two values Am 2 ol = O(10~ 5 eV 2 ) indicated in the first well as 

the value Am 2 ol = O(10 -10 eV 2 ) allowed in the second are situated — consistently - 
much below the region excluded for disappearance modes of v e by CHOOZ experiment, 
where Ara 2 CH ~ 0.9 x 10~ 3 eV 2 (at the 90% confidence level) [6]. 

We can see from the experimental solar-neutrino estimate (49) that it may be related 
to our neutrino oscillation formula (28) (both being likely enhanced by the MSW mech- 
anism). Then, only the disappearance mode v e — > is s , described by the formula (27), 
contributes to the rhs of Eq. (20). In such a case we get with the use of figures (50) two 
other neutrino inputs 



and 



4F 2 

1 + y2)2 = sin2 2d «* ~ 8 x 10 - 3 ( 53 ) 



|m 2 4 -m 2 J = Am s 2 ol ~ 5 x 10~ 6 eV 2 , (54) 
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if the first MSW solution is taken into account. The input (53) gives Y ~ 0.045 (and 
Y 2 ~ 2.0 x 10~ 3 ). 

From the first Eq. (10) with M$ = we obtain 



I M {v) I 

Y ~ (55) 

under the assumption that M^ 2 > (2\M$\) 2 . Thus, (\M$ l/M^f ~ 0.0020 with the 
input (53). 

On the other hand, from the first Eq. (8) with M$ = 

-1W44 ■'"44 

when Afff 2 > (2|M 1 ( 4 ) |) 2 . Thus, M$ 2 ~5x 10~ 6 eV 2 with the input (54), and so, 
\M$\ 2 ~ 1.0 x 10" 8 eV 2 due to the input (53). Hence, 

m vi 4.5 x 10~ 6 eV , m Vi ~ 2.2 x 10~ 3 eV . (57) 

Here again, the minus sign at m Ul is phenomenologically irrelevant. 

We can see from Eq. (20) that its approximate form (28) valid for (a^/m^ — > 0, 
applied here, works very well even for the central value of (a^/m^ , because then 

4F 2 „ _ , , 16 (a^ 2 



8 x 10" 3 large versus ~ 2.5 x 10~ 4 . (58) 

(1 + F 2 ) 2 & 841 \mj K ' 

In fact, due to Eq. (5), 

16 fa^\ 2 

< ^ _ < 6.2 x 10~ 4 , (59) 
" 841 \mj ~ K J 

where the central value is 2.5 x 10" 4 . 

The small value (59) is the main reason, why in the case of no sterile neutrino v s 

our neutrino oscillation formula (20) cannot be compared with the experimental estimate 

(49). Indeed, in this case it assumes the form 

P („ e - „ e ) = 1 - Jl (sin 2 x 21 + X 2 sin 2 s 3 i) • (60) 
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To conclude the above discussion of solar neutrinos, we can claim that, in the frame- 
work of our "texture", only the disappearance mode of u e to the sterile neutrino v s (likely 
to be enhanced by the MSW mechanism in the Sun matter) can be responsible for the 
observed deficit of solar z/ e 's. The predicted neutrino masses are \m Ul \ ~2x 1CT 6 eV and 
m VA ~ 2 x 10~ 3 eV with m 2 4 - m 2 vi ~ 5 x 1(T 6 eV 2 . 

Our last remark concerns the LSND experiment that seems to detect — > u e oscil- 
lations by observing the appearance of v e originating from produced in n + decay [9]. 
The observed excess of z/ e 's, analized in terms of two-family oscillation formula, leads 
to sin 2 26 l L sND > 4 x 10" 4 , in particular to sin 2 2# LSND ~ 1.5 x 10~ 3 - 1.5 x KT 1 for 
Am1 SND ~ 1 eV 2 (c/. [9]). This shows that, due to the small value (59), the LSND mag- 
nitude of — > v e oscillation amplitude can hardly be explained by means of our formula 
(17) for P(Vfj, — > v e ) = P(y e — > v^) . However, there is possibly a narrow overlap near 
the upper limit of the range (59), corresponding to the values sin 2 26> L snd ~ 6 x 10~ 4 and 
Am 2 SND ~ 1.2 eV 2 , located at the border of LSND allowed region (at the 95% confidence 
level) [9]. Unfortunately, another (earlier) LSND experiment on — > v e oscillations es- 
timates that an analogical allowed region is much narrower (c/. [9]), what excludes the 
above potential overlap. 

Nevertheless, it may be interesting to remark that, in the case of values m 2 2 ~ m 2 3 ~ 
1.2 eV 2 corresponding to the large Am 2 SND ~ 1.2 eV 2 , we can put xi\ ~ ^31 ~ ^24 ~ 
0:34 3> a; 32 > £41 in Eq. (17). This leads to the two-family oscillation formula 



when short-baseline experiments with x 2 i ~ tt/2 are considered (here, 1 + Y 2 ~ 1 + 
0.0020 ~ 1). 

6. Outlook: an unconventional picture 

If beside the sterile neutrino v s there exists also a second sterile neutrino mentioned in 
the footnote * , two opposite options seem to be attractive. The first was already outlined 
in this footnote: such a second sterile neutrino might mix mainly with leading to a new 
disappearance mode of that, if not too strong, would only slightly correct the effect of 
the dominating mode — > u T . In the second option, the roles of these two modes of 




(61) 
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would be interchanged: the new disappearance mode of would dominate the mode — > 
u T1 and so would be responsible for producing a near-to-maximal oscillation amplitude 
for — > as is observed in Super-Kamiokande. In such a case, both parameters 
and (3^ might be small, perhaps zero. 

Thus, this second option would create a uniform but unconventional picture of neutrino 
oscillations, where they would be caused essentially by mixing two sterile neutrinos with v e 
and z/^, respectively. At any rate, this picture would be true for the solar and atmospheric 
neutrinos. 

Sterile neutrinos of both kinds might constitute an important part of relativistic dark 
matter, passive with respect to all Standard-Model interactions (including weak interac- 
tions). In such a case, they might even be the main constituents of matter in the Universe. 
Of course, sterile neutrinos would interact with each other and with Standard-Model ac- 
tive particles through gravitation. 

Appendix 

For the neutrino mass matrix (M^) as given in Eq. (1) we assumed thatch = [Eq. 
(2)], what implied M$ = =M^ ] . As we saw, this assumption, leading to weak mixing 
oiu e withz/^ andz/ T , is neatly consistent with negative results of CHOOZ experiment [6]. 
Now, we will relax such an extremal assumption, allowing for a nonzero but small a^, 
much smaller than the large (3^ = O(10 -4 to 1) eV [Eq. (48)] responsible for near- 
to-maximal mixing of with u T , just as suggested by results of Super-Kamiokande 
experiment [4,5]. Although, in our neutrino "texture" such a relaxation is not needed to 
understand solar neutrino experiments which are here reasonably explained by the mixing 
of v e with v s , it may be applied to appearance experiments in the mode — > v e or v e — > v^. 

In this case, (m\^ is perturbed by the matrix 



K) 



29 
Thus, 



( 2e^ M \ 

2e-^ M SVSe' 1 ^ 

Zyfar**™ j 



(A.l) 
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5M\ v) 



a 



M 



"12 



29 



<5M 9 (iy)! 



'21 



<5M 9 H 



'23 



29 



'32 



while 



(A.2) 







M 2 <? 



M. 



M 



23 



29 ' 



M. 



M> 



32 



(A.3) 



Then, the total secular equation det + SM^ — l(m Vi + 5m v% ) = gives in the 
lowest (quadratic or linear) perturbative order in / fi^ the following neutrino mass 
corrections: 



<5m 



5m 



"2, V3 



Ml 2 jitW 



<5M^| 2 M 



33 



'12 



1 - 



|M<?| 



;/ 2 



^3 



2|5M 2 ( 3 ) ||M 



= ± 



rMi 

23 I 



X 



771 



m 



V, 



= ±2|5Mg ) | IT ^ ) 



X 



(A.4) 



where \6M$\ = 0.0690c*M \6M$\ = 0.478aM, M$ = 20.7/i^ and \M$\ = 0.478 
When the mass matrix (m,-^) is perturbed by the matrix ((5M^) given in Eq. (A.l), 



then the unitary diagonalizing matrix (U^\ [Eq.(9)] undergoes the perturbation [5U l 



r(") 



which, after some calculations in the lowest perturbative order, can be written as follows: 



SU. 



M 





\6Mg>\M$ c _ iv M 



\sm£>\ 



V 



with m!# /m U3 = 1 — \M^\X/m Vs and m V2 = —\m U2 \. Here, Y ~ 0.0045 is put zero for 
simplicity, and i , j — 1, 2, 3 (if Y ^ 0, then 511^ and ^[/^ get extra factors 1/VT+ F 2 ). 
In particular, the elements (Jt/^ an d ^31^ are produced (in the lowest perturbative 
order) through the diagonalizing procedure of the total mass matrix (m^ + 5M-^ in 
the following way: 



Mi 



2i V <"' 





|MM|(l+X2)3/2 



|5A^>|X(1-X 2 ) 
|M^|(l+X2)3/2 







(A.5) 



/ 



utf + su. 



M 
21 



f/ 3 ( l ) +^: 



M 

31 



m, 



<5m 



^1 



»/fM 1 t\A v ) 
\5M 12 \M 33 __ i(p („) 



m,. 



(Jm^ ) (M^ + 6M&> ) <Jm^ 



(M& } + <5M^ ) (M^ - m^ - <Jm^ ) <5M^ M ; 



r M> 



r(") 



M *A») 



33 



771 U2 777 



23 I -2i ¥ )(") 

16 ' 



(A.6) 



where Eq. (A. 4) for 5m Ul is used, while f/^i — an d ^31 — 0- Of course, M^Y = and 
so, m vi = for Y = 0. 

Once (SU-j^ is known, the perturbation (SVij) of the leptonic Cabibbo — Kobayashi— 
Maskawa matrix (V^-) [Eqs. (15)] can be evaluated from the definition: 

Pa + SVy) = (t#> + (Ujg) . (A.7) 

Here, u a = J2j (v* a + 5V*Jj vj with v a = v e , , v T and vj = i>i, v 2 , v z . Thus, from Eq. 
(A.7) we obtain 

W 3 ) = (U<$) * , (A.8) 

where, in the second step, terms proportional to {a^ /m U2 ){a^ /m^) are neglected versus 
terms proportional to (aM/m^). Thus, in this case 



SVij = SU*, . (A.9) 

Making use of the elements SV^ defined through Eqs. (A.9) and (A. 5), we can calculate 
from Eqs. (16) the perturbations 5P{v a — > i/g) of the neutrino-oscillation probabilities 
P{y a — > 1/3) which were evaluated before in the case of M$2 = = M 2 ( i } [Eqs. (17) - 
(22)]. 

In particular for u e — > 1^ oscillations, after some calculations up to the quadratic 
perturbative order in aW/juW, we obtain the following formula correcting Eq. (17): 



8P(v e -> z/ M ) ~ 



16 aMa( e ) [ M 3 ( 3 } 



841 I I I 



sin 2 (a;2i + <p) + sin 2 (0:21 — v?) — 2 sin 2 (p 
X 2 



16 /aM' 



(1 + X 2 ) 5 
X 2 



sin (x 3 2 + v?) - sin (x 32 - V?) 



841 UJ (TTW sm X32 ' 



(A.10) 
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where cp = (1/2) — ip^). Here, Y = for simplicity, and the upper limit of X ~ 0.999 
is taken into account. Owing to this, in our calculations leading to Eq. (A. 10), the term 
proportional to {M<$ /m V3 )(\ — X 2 \m 1/2 \/m„ 3 ) ~ 9.04 x 10~ 6 is negligible and the relation 
X21 ~ X31 ^ ^32 works. In fact, for such a value of X, we have \m U2 \ ~ 1.05 eV ~ m V3 
(m U2 = — \m V2 \) with m 2 3 — m 2 2 ~ 5 x 10~ 3 eV 2 and \m V2 \/m vz ~ 0.998 ~ 1, while 
~ 1.08 x 10~ 4 eV and /jM ~ 2.20 eV with (3^/m u . A ~ 2.09. Further, 



^ 2 , 2 x 10 -3 „ * l^M. , _J_ „ . 5Q1 , „ . 499 (A . n) 

and, from Eqs. (A. 4), 

, /«W\ 2 «(") 

tfm^ ~ -1.06 x 10" 5 eV , 5m U2 „ 3 ~ ±0.504 eV . (A.12) 

\m U2 J ' 

Since the second term in Eq. (A. 10) is not invariant under the change of phase sign, 
ip — > —ip, this term violates time reversal and so, CP reflection (as CPT is conserved). 
The formula (A. 10) can be rewritten in the following numerical form: 



M 33 



5P(v e -> z/ M ) ~ [4.65 x 10~ 6 [sin 2 (a;2i + <p) + sin 2 (a;2i -</?)- 2 sin 2 y9 



- 5.47 x 10~ 4 



sin 2 (x 32 + ip) - sin 2 (a;32 - <p)] } 
+4.76 x 10~ 3 ( — ) sin 2 x 32 • (A. 13) 



If ip — > and a^/l 772 ^! < 1> the second and third term here can be neglected, when 
short-baseline experiments with x 2 i ~ 7r/2 are taken into account, since then sin 2 x 32 — 
x\ 2 ~ 2.5 x 10- 5 (tt/2) 2 . In such a case, therefore, we get the two-family formula 

5P(z/ e -> ~ 9.30 x 10 -6 - sin 2 x 21 . (A.14) 

\TfLy 2 I 

Evidently, the oscillation amplitudes in the perturbative formula (A. 13) are much too 
small to be able to help the unperturbed formula (61) in explaining the LSND estimate 
for — > z/ e oscillations [9]. 
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